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The first reports of and ?/>' production from experiment 
NA50 at the CERN SPS are compared to calculations based 
on a hadronic model of charmonium suppression developed 
previously. Data on centrality dependence and total cross 
sections are in good accord with these predictions. 



Experiment NA50 has reported an abrupt decrease in 
ib production in Pb+Pb collisions at 158 GeV per nucleon 
jjj. Specifically, the collaboration presented a striking 
'threshold effect' in the i/j-to-continuum ratio by plot- 
ting it as a function of a calculated quantity, the mean 
path length of the ib through the nuclear medium, L, as 
shown in fig. la. This apparent threshold has sparked 
considerable excitement as it may signal the formation 
of quark-gluon plasma in the heavy Pb+Pb system . 
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FIG. 1. (a) The NA50 § comparison of ib production in 
Pb+Pb and S+U collisions as a function of the average path 
length L, see eq. (3). B is the ip — > p + fi~ branching ratio, 
(b) Transverse energy dependence of Pb+Pb data. Curves in 
(a) and (b) are computed using eqs. (4-6). 

Our aim is to study the Pb results in the context of 
a hadronic model of charmonium suppression HQ. We 
first demonstrate that the behavior in the NA50 plot, 
fig. la, is not a threshold effect but, rather, reflects the 
approach to the geometrical limit of L as the collisions 
become increasingly central. When plotted as a func- 
tion of the measured neutral transverse energy Et as in 
fig. lb, the data varies smoothly as in S+U measurements 
in fig. 3b below The difference between S+U and 

Pb+Pb data lies strictly in the relative magnitude. To 
assess this magnitude, we compare if) and ib' data to ex- 



pectations based on the hadronic comover model §,§. 
The curves in fig. 1 represent our calculations using pa- 
rameters fixed earlier in Ref. Q . Our result is essentially 
the same as the Pb+Pb prediction in || . 

The hadronic contribution to the suppression arises 
from scattering of the nascent ib with produced parti- 
cles - the comovers - and nucleons HQ. To determine 
the suppression from nucleon absorption of the ib, we cal- 
culate the probability that a cc pair produced at a point 
(6, z) in a nucleus survives scattering with nucleons to 
form a ib. The standard MM result is 



(1) 



S A = cxp{- / dzp A (b, z)(t^ n } 

J z 



where pa is the nuclear density, b the impact parame- 
ter and o~ip]\[ the absorption cross section for ?/;-nucleon 
interactions. One can estimate S A ~ exp{— ct^atPo^a}, 
where L A is the path length traversed by the cc pair. 

Suppression can also be caused by scattering with 
mesons that happen to travel along with the cc pair (see 
refs. in ||). The density of such comovers scales roughly 
as Et- The corresponding survival probability is 



5 co = exp{- / dTna co v re i} 



(2) 



where n is the comover density and r is the time in the 
ib rest frame. We write S co ~ exp{— f3Er}, where /? 
depends on the scattering frequency, the formation time 
of the comovers and the transverse size of the central 
region, R T , cf. eq. (8). 

To understand the saturation of the Pb data with L in 
fig. la, we apply the schematic approximation of Ref. 
for the moment to write 
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where the brackets imply an average over the collision 
geometry for fixed Et and <j(Et) = da/dEr- The path 
length L = (La + Lb) and transverse size Rt depend on 
the collision geometry. The path length grows with Et, 
asymptotically approaching the geometric limit Ra+Rb- 
Explicit calculations show that nucleon absorption begins 
to saturate for b < Ra, where Ra is the smaller of the two 
nuclei, see fig. 4 below. On the other hand, Et continues 
to grow for b < Ra due, e.g., to fluctuations in the num- 
ber of NN collisions. Equation (2) falls exponentially in 
this regime because (3, like L, saturates. 



1 



In fig. lb, we compare the Pb data to calculations of 
the ^-to-continuum ratio that incorporate nucleon and 
comover scattering. The contribution due to nucleon ab- 
sorption indeed levels off for small values of b, as ex- 
pected from eq. (3) . Comover scattering accounts for the 
remaining suppression. 

These results are predictions obtained using the com- 
puter code of Ref. J3| with parameters determined in 
Ref. 0]. However, to confront the present NA50 analysis 
0, we account for changes in the experimental coverege 
as follows: 



We parametrize P(E T ,b) = Cexp{-(E T - E T ) 2 /2A}, 
where E T (b) = eA/"(&), A(b) = ueE T (b), C(b) = 
(27rA(6)) _1 and Af(b) is the number of participants (see, 
e.g., Ref. ||). We take e and to to be phenomenological 
calorimeter-dependent constants. 

We compare the minimum bias distributions for total 
hadronic Et calculated using eq. (7) for e — 1.3 GeV and 
lu = 2.0 to NA35 S+S and NA49 Pb+Pb data @. The 
agreement in fig. 2a builds our confidence that eq. (7) 
applies to the heavy Pb+Pb system. 



• Calculate the continuum dimuon yield in the new 
mass range 2.9 < M < 4.5 GeV. 

• Adjust the Et scale to the pseudorapidity accep- 
tance of the NA50 calorimeter, 1.1 < rj < 2.3. 

The agreement in fig. 1 depends on these updates. 

We now review the details of our calculations, high- 
lighting the adjustments as we go. For collisions at a 
fixed b, the ^-production cross section is 
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(b) 



a$ N I d 2 sdzdz' p A (s, z)p B {b - s, z') S, (4) 



where S = SaSbS co is the product of the survival proba- 
bilities in the projectile A, target B and comover matter. 
The continuum cross section is 
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The magnitude of (4,5) and their ratio are fixed by the 
elementary cross sections a^ N and 0^+^- ■ We calculate 
o~^ N using the phenomenologically-successful color evap- 
oration model p0|| . The continuum in the mass range 
used by NA50, 2.9 < M < 4.5 GeV, is described by 
the Drell-Yan process. To confront NA50 and NA38 
data in the appropriate kinematic regime, we compute 
these cross sections at leading order following jlO 11 1 us- 
ing GRV LO parton distributions with a charm if-factor 
K c = 2.7 and a color evaporation coefficient — 2.54% 
and a Drell-Yan if-factor Kdy — 2.4. Observe that 
these choices were fixed by fitting pp data at all available 
energies fo). Computing a£f for 2.9 < M < 4.5 GeV 
corresponds to the first update. 

To obtain Et dependent cross sections from eqs. (4) 
and (5), we write 



AB 



(E T ) = Jd 2 bP(E T ,b)a AB 



(b). 



(6) 



The probability P(ET,b) that a collision at impact pa- 
rameter b produces transverse energy Et is related to the 
minimum-bias distribution by 



.(Et) = J d 2 b P(E 7 



,b). 



(7) 
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FIG. 2. Transverse energy distributions from eq. (7). The 
S-Pb comparison (a) employs the same parameters. 

Figure 2b shows the distribution of neutral transverse 
energy calculated using eqs. (5) and (6) to simulate the 
NA50 dimuon trigger. We take e = 0.35 GeV, u) = 3.2, 



NN 



37.2 pb as appropriate for the dimuon- 



mass range 2.9 < M < 4.5 GeV. The E T distribu- 
tion for S+U — > p + p~ + X from NA38 was described 
using e = 0.64 GeV and u> = 3.2 - the change in 
e corresponds roughly to the shift in particle produc- 
tion when the pseudorapidity coverage is changed from 
1.7 < 7) < 4.1 (NA38) to 1.1 < 77 < 2.3 (NA50). Tak- 
ing e = 0.35 GeV for the NA50 acceptance is the second 
update listed earlier. 

We now apply eqs. (1,2,4) and (5) to charmonium sup- 
pression in Pb+Pb collisions. To determine nucleon ab- 
sorption, we used pA data to fix a^N ~ 4.8 mb in Ref. [Q. 
This choice is in accord with the latest NA38 and NA51 
pA data, see fig. 3a. To specify comover scattering [Q, 
we assumed that the dominant contribution to tp dissoci- 
ation comes from exothermic hadronic reactions such as 
p + ip — > D + D. We further took the comovers to evolve 
from a formation time tq ~ 2 fm to a freezeout time 
Tp ~ Rt/v 1c i following Bjorken scaling, where v vc i ~ 0.6 
is roughly the average i/j— p relative velocity. The survival 
probability, eq. (2), is then 



S co = CXp{-O- co ti re l7l0T'0 ln(i?T/VrclTo)} 



(8) 



where a co ~ 2er^j\r/3, Rt ~ Ra and no is the initial 
density of sufficiently massive p, lu and r\ mesons. To 
account for the variation of density with Et, we take 
no = tioEt/Et(0) B- A value n — 0.8 fm~ 3 was chosen 
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to fit the central S+U datum. Since we fix the density 
in central collisions, this simple ansatz for S co may be 
inaccurate for peripheral collisions. [Densities ~ 1 fm 
typically arise in hadronic models of ion collisions, e.g., 
rcfs. [jl4j. The internal consistency of hadronic models at 
such densities demands further study] 

We expect the comover contribution to the suppres- 
sion to increase in Pb+Pb relative to S+U for central 
collisions because both the initial density and lifetime of 
the system can increase. To be conservative, we assumed 
that Pb and S beams achieve the same mean initial den- 
sity. Even so, the lifetime of the system essentially dou- 
bles in Pb+Pb because Rt ~ Ra increases to 6.6 fm 
from 3.6 fm in S+U. The increase in the comover contri- 
bution evident in comparing figs, lb and 3b is described 
by the seemingly innocuous logarithm in eq. (8), which 
increases by w 60% in the larger Pb system. 



would not appreciably change our result. Note that a 
uniform decrease of the ratio would increase the comover 
contribution needed to explain S+U collisions. 

To see why saturation occurs in Pb+Pb collisions but 
not in S+U, we compare the NA50 L(Et) B to the av- 
erage impact parameter (6}(£t) in fig- 4. To best un- 
derstand fig. la, we show the values of L(Et) computed 
by NA50 for this figure. We use our model to compute 
(6) = (bT AB )/(T AB ), where (/(&)) = Jd 2 b P(E T ,b)f(b) 
and Tab = Jdrsdzdz' p A (s, z)ps{b — s,z'). [Note that 
NA50 reports similar values of (6) (St) PJ.] In the 
Et range covered by the S experiments, we see that 
(b) is near ~ i?s = 3.6 fm or larger. In this range, 
increasing b dramatically reduces the collision volume 
and, consequently, L. In contrast, in Pb+Pb collisions 
(6) <C Rpb = 6.6 fm for all but the lowest Et bin, so that 
L does not vary appreciably. 
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FIG. 3. (a) pA cross sections [^j] in the NA50 acceptance 
and (b) S+U ratios from '91 § and '92 jjj runs. The '92 data 
are scaled to the '91 continuum. The dashed line indicates the 
suppression from nucleons alone. The pp cross section in (a) 
is constrained by the global fit to pp data in ref. jl^ . 

In Ref. [El, we pointed out that comovers were nec- 
essary to explain S+U data from the NA38 1991 run 
||. Data just released [jO from their 1992 run support 
this conclusion. The '91 ip data were presented as a 
ratio to the dimuon continuum in the low mass range 
1.7 < M < 2.7 GeV, where charm decays are an im- 
portant source of dileptons. On the other hand, the '92 
ip data jl]]8| are given as ratios to the Drell-Yan cross 
section in the range 1.5 < M < 5.0 GeV. That cross 
section is extracted from the continuum by fixing the K- 
factor in the high mass region ju| . To compare our result 
from Ref. Q to these data, we scale the '92 data by an 
empirical factor. This factor is w 10% larger than our 
calculated factor a%$(92)/a%£ t (91) « 0.4; these values 
agree within the NA38 systematic errors. [NA50 simi- 
larly scaled the '92 data to the high-mass continuum to 
produce fig. la.] Because our fit is driven by the highest 
Et datum, we see from fig. 3b that a fit to the '92 data 
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FIG. 4. E T dependence of L (solid) used by NA50 § (see 
fig. la) and the average impact parameter (6) (dot-dashed). 
The solid line covers the measured Et range. 

NA50 and NA38 have also measured the total ip- 
production cross section in Pb+Pb Q and S+U reactions 
||. To compare to that data, we integrate eqs. (4, 6) to 
obtain the total (a/AB)^ = 0.95 nb in S+U at 200 GeV 
and 0.54 nb for Pb+Pb at 158 GeV in the NA50 spec- 
trometer acceptance, 0.4 > xjr > and —0.5 < cos# < 
0.5 (to correct to the full angular range and 1 > xp > 0, 
multiply these cross sections by « 2.07). The experi- 
mental results in this range are 1.03 ± 0.04 ± 0.10 nb for 
S+U collisions § and 0.44 + 0.005 + 0.032 nb for Pb+Pb 
reactions [Q. Interestingly, in the Pb system we find a 
Drell-Yan cross section {a/AB) DY = 37.2 pb while NA50 
finds (a / AB) DY = 32.8 ± 0.9 ± 2.3 pb. Both the V and 
Drell-Yan cross sections in Pb+Pb collisions are some- 
what above the data, suggesting that the calculated rates 
at the NN level may be - 20-30% too large at 158 GeV. 
This discrepancy is within ambiguities in current pp data 
near that low energy fl0]| . Moreover, nuclear effects on 
the parton densities omitted in eqs. (4,5) can affect the 
total S and Pb cross sections at this level. 
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We remark that if one were to neglect comovers and 
take (T^at = 6.2 mb, one would find (a/AB)^ = 1.03 nb 
in S+U at 200 GeV and 0.62 nb for Pb+Pb at 158 GeV. 
The agreement with S+U data is possible because co- 
movers only contribute to the total cross section at the 
~ 18% level in the light system. This is expected, since 
the impact-parameter integrated cross section is domi- 
nated by large b and the distinction between central and 
peripheral interactions is more striking for the asymmet- 
ric S+U system. As in Ref. 0, the need for comovers is 
evident for the £/r - dependent ratios, where central col- 
lisions are singled out. 
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show that the singlet and octet extrapolations describe 
pA data equally well. 

Our predictions for Pb+Pb collisions are shown in 
fig. 6. In the octet model, the entire suppression of the 
ip'-to-ip ratio is due to comover interactions. In view of 
the schematic nature of our approximation to S co in eq. 
(8), we regard the agreement with data of singlet and 
octet extrapolations as equivalent. 

In summary, the Pb data |lj cannot be described by 
nucleon absorption alone. This is seen in the NA50 plot, 
fig. la, and confirmed by our results. The saturation with 
L but not Ex suggests an additional density-dependent 
suppression mechanism. Earlier studies pointed out that 
additional suppression was already needed to describe 
the S+U results [Q; recent data 0] support that con- 
clusion (see, however, Q). Comover scattering explains 
the additional suppression. Nevertheless, it is unlikely 
that this explanation is unique. SPS inverse-kinematics 
experiments (B < A) and AGS pA studies near the ip 
threshold can help pin down model uncertainties. 

This work was supported in part by US-DOE contracts 
DE-AC02-76CH00016 and DE-AC03-76SF0098. 

Note Added — After the completion of this manuscript, 
we learned of cascade calculations |Q that confirm our 
conclusions. Such calculations do not employ the simpli- 
fications (e.g. no oc Et) needed to derive (8). Some of 
these authors took ct^jv ~ 6 mb (instead of ~ 5 mb) to 
fit the NA51 data in fig. 3a somewhat better. 
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To apply eqs. (4-6) to calculate the ip'-to-ip ratio as a 
function of Et, we must specify o~^, N , o~^>n, and <j^i co . 
Following Ref. ffij , we use pp data to fix Bo-$, N /Ba$ N = 
0.02 (this determines F^i). The value of ct^n depends 
on whether the nascent ip 1 is a color singlet hadron or 
color octet cc as it traverses the nucleus. In the singlet 
case, one expects the absorption cross sections to scale 
with the square of the charmonium radius. Taking this 
ansatz and assuming that the if)' forms directly while 
radiative x decays account for 40% of ip production, one 
expects o~ipt ~ 2.1cr^, for interactions with either nuclcons 
or comovers For the octet case, we take cr^'N ~ o~^n 
and fix o~f> co w 12 mb to fit the S+U data. In fig. 5a, we 
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